Teaching the Kepler laws for freshmen. 
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Abstract 

We present a natural proof of Kepler's law of ellipses in the spirit 
of Euclidean geometry. Moreover we discuss two existing Euclidean 
geometric proofs, one by Feynman in hist Lost Lecture from 1964 and 
the other by Newton in the Principia of 1687. 

1 Introduction. 

One of the highlights of classical mechanics is the mathematical derivation 
of the three experimentally observed Kepler laws of planetary motion from 
Newton's laws of motion and of gravitation. Newton published his theory 
of gravitation in 1687 in the Principia Mathematica [8]. After two short 
introductions, one with definitions and the other with axioms (the laws of 
motion), Newton discussed the Kepler laws in the first three sections of Book 
1 (in just 40 pages, without ever mentioning the name of Kepler!) 

Kepler's second law (motion is planar and equal areas are swept out in 
equal times) is an easy consequence of the conservation of angular momentum 
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L = r X p, and holds in greater generality for any central force field. All this 
is explained well by Newton in Propositions 1 and 2. 

Kepler's first law (planetary orbits are ellipses with the center of the force 
field at a focus) however is peculiar for the attractive force field. Using 
Euclidean geometry Newton derives in Proposition 11 that the Kepler laws 
can only hold for an attractive 1/r^ force field. The reverse statement that 
an attractive 1/r^ force field leads to eUiptical orbits Newton concludes in 
Corollary 1 of Proposition 13. Tacitly he assumes for this argument that the 
equation of motion F = ma has a unique solution for given initial position 
and initial velocity. Theorems about existence and uniqueness of solutions of 
such a differential equation have only been formulated and mathematically 
rigorously been proven in the 19th century. However there can be little doubt 
that Newton did grasp these properties of his equation F = ma [1]. 

Somewhat later in 1710 Jakob Hermann and Johan Bernoulli gave a direct 
proof of Kepler's first law, which is still the standard proof for modern text 
books on classical mechanics |9]. One writes the position vector r in the 
plane of motion in polar coordinates r and 9. The trick is to transform the 
equation of motion ma = —kr/r^ with variable the time t into a second order 
differential equation of the scalar function u = 1/r with variable the angle 6. 
This differential equation can be exactly solved, and yields the equation of 
an ellipse in polar coordinates [3]. 

Another popular proof goes by writing down the so called Laplace-Runge- 
Lenz vector 

K = p X L — kmr/r 

with p = mv momentum and F(r) = —kr/r^ the force field of the Kepler 
problem. The LRL vector K turns out to be conserved, i.e. K = 0. This 
result can be derived by a direct computation as indicated in Section 2. An 
alternative geometric argument is sketched in Section 3. Working out the 
equation r ■ K = rK cos 6 gives the equation of an ellipse in polar coordinates 
[3]. The geometric meaning of the LRL vector becomes clear only a posteriori 
as a vector pointing in the direction of the long axis of the ellipse. But the 
start of the proof to write down the LRL vector remains a trick. For a 
historical account of the LRL vector we refer to Goldstein [1]. 

The purpose of this note is to present in Section 2 a proof of the Kepler 
laws for which a priori the reasoning is well motivated in both physical and 
geometric terms. In Section 3 we review the hodographic proof as given by 
Feynman in his "Lost Lecture" [5], and in Section 4 we discuss Newton's 
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proof from the Principia[H]. All three proofs are based on Euclidean geom- 
etry, although we do use the language of vector calculus in order to make 
it more readable for people of the 21st century. We feel that our proof is 
really the simplest of the three, and at the same time it gives more refined 
information (namely the value of the long axis 2a = —k/H of the ellips 
£). In fact we think that our proof in Section 2 can compete equally well 
both in transparency and in level of computation with the standard proof of 
Jakob Hermann and Johann Bernoulli, making it an appropriate alternative 
to present in a freshman course on classical mechanics. 

We like to thank Hans Duistermaat and Arnoud van Rooij for useful 
comments on this article. 



2 A Euclidean proof of Kepler's first law. 

We shall use inner products u ■ v and outer products u x v of vectors u and 
V in with the compatibility conditions 

u • (v X w) = (u X v) ■ w 

u X (v X w) = (u ■ w)v — (u ■ v)w 
and the Leibniz product rules 

(u ■ v)' = U ■ V + U ■ V 

(u X v)' = U X V + U X V 

without further explanation. 

For a central force field F(r) = /(r)r/r the angular momentum vector 
L = r X p is conserved by Newton's law of motion F = p, thereby leading to 
Kepler's second law. For a spherically symmetric central force field F(r) = 
/(r)r/r the energy 

H = p'^/2m + V{r) , V{r) = - j f{r) dr 

is conserved as well. These are the general initial remarks. 

From now on consider the Kepler problem /(r) = —kjr^ en Vir) = —k/r 
with > a coupling constant. For fixed energy H < the motion is 
bounded inside a sphere with center and radius —k/H. Consider the 
following picture of the plane perpendicular to L. 
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The circle C with center and radius —k/H is the boundary of a disc 
where motion with energy H < takes place. Let s = —kr/rH be the 
projection of r from the center on this circle C. The line C through r with 
direction vector p is the tangent line of the orbit 8 at position r with velocity 
V. Let t be the orthogonal reflection of the point s in the line £. 

Theorem. The point t equals K./mH and therefore is conserved. 

Proof. The line J\f spanned by n = p x L is perpendicular to C The point 
t is obtained from s by subtracting twice the orthogonal projection of s — r 
on the line J\f, and therefore 

t = s - 2((s - r) • n)n/n^ 

Now 

s = —kr/rH 
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(s - r) ■ n = -(if + k/r)T ■{pxL)/H = -{H + k/r)L'^/H 
= p^L^ = 2m{H + k/r)L^ 

and therefore 

t = —kr/rH + n/mH = K./mH 

with K = p X L — kmr/r the LRL vector. The final step K = is derived 
by a straightforward computation using the compatibihty relations and the 
Leibniz product rules for inner and outer products of vectors in R'^. □ 

Corollary. The orbit E is an ellipse with foci and t, and long axis equal 
to 2a = -k/H. 

Proof. Indeed we have 

|t - r| + |r - 0| = |s - r| + |r - 0| = |s - 0| = -k/H. 

Hence S is an ellipse with foci and t, and long axis 2a = —k/H . □ 

The above proof has two advantages over the earlier mentioned proofs of 
Kepler's first law. The conserved vector t = K/mif is a priori well motivated 
in geometric terms. Moreover we use the gardener's definition of an ellipse. 
The equation of an ellipse in polar coordinates is unknown to most freshmen, 
and so additional explanation would be needed for that. Another advantage 
of our proof is that the solution of the equation of motion is achieved by just 
finding enough constants of motion. The proofs by Feynman and Newton in 
the next sections on the contrary rely at a crucial point on the existence and 
uniqueness theorem of differential equations. 

We proceed to derive Kepler's third law along standard lines [3]. The 
ellipse 8 has numerical parameter^ a, 6, c > given by 2a = —k/H, 4c^ = 
K^/rn^H^ = {2mHL'^ + k"^) / rn^ H"^ and = 6^ + c^. The area of the 
region bounded by S equals 

TTab = LT /2m 

with T the period of the orbit. Indeed L/2m is the area of the sector swept 
out by the position vector r per unit time. A straightforward calculation 
gives 

TVa^ = An^m/k. 
^The long axis equals 2a, the short axis 2b, and = + c^. 
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Newton's law of gravitation states that the couphng constant k is propor- 
tional to the product of the mass m of the planet and the mass M of the sun, 
and we find Kepler's third (harmonic) law stating that T"^ /a? is the same for 
all planet^. 

3 Feynman's proof of Kepler's first law. 

In this section we discuss a different geometric proof of Kepler's first law 
based on the hodograph Ti. By definition Ti is the curve traced out by the 
velocity vector v in the Kepler problem. This proof goes back to Mobius 
in 1843 and Hamilton in 1845 [2j and has been forgotten and rediscovered 
several times, among others by Feynman in 1964 in his "Lost Lecture" [S]. 

Let us assume (as in the picture of the previous section) that ivn/n = v 
with i the counterclockwise rotation around over tt/2. So the orbit S is 
assumed to be traversed counterclockwise around the origin 0. 

Theorem. The hodograph Ti is a circle with center c = i'K./mL and radius 
k/L 

Proof. We shall indicate two proofs of this theorem. The first proof is analytic 
in nature, and uses conservation of the LRL vector K by rewriting 

K = p X L — kmr/r = mvLn/n — kmv/r 

as 

vn/n = K/mL + kv/rL , 

or equivalently 

V = iK./mL + ikr/rL . 

Hence the theorem follows from K = 0. 

There is a different geometric proof of the theorem, which as a corollary 
gives the conservation of the LRL vector K. The key point is to reparametrize 
the velocity vector v from time t to angle 6 of the position vector r. It turns 

■^The mass m we have used so far is in fact equal to the reduced mass mM/(rn + M), 
and this almost equals m if m -C Af. The coupling constant k is according to Newton 
equal to GmM with G the universal gravitational constant. We therefore see that Kepler's 
third law holds only approximately for m <^ M . 
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out that the vector v(^) is traversing the hodograph H. with constant speed 
k/L. Indeed we have from Newton's laws 



(iv dt kr dt 

and Kepler's second law gives 



r^d9/2 = Ldt/2m . 
Combining these identities yields 

so indeed v(^) travels along H with constant speed k/L. Since r = re*^ a 
direct integration yields 

v(^) = c + ikr/rL , c — 

and the hodograph becomes a circle with center c and radius k/L. Compar- 
ison with the last formula in the first proof gives 

c = iK/mL 

and K = comes out as a corollary. □ 

All in all, the circular nature of the hodograph H is more or less equivalent 
to the conservation of the LRL vector K. 
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Now turn the hodograph Ti. clockwise around over tt/2 and translate 
over ic = — K/mL. This gives a circle V with center and radius k/L. Since 

kr/rL + K/mL = vn/n = — iv 

the orbit S intersects the line through and kr/rL in a point with tangent 
line C perpendicular to the line through kr/rL and — K/mL. For example 
the ellipse JF with foci and —K/mL and long axis equal to k/L has this 
property, but any scalar multiple AJF with A > equally does so. Since 
curves with the above property are uniquely charcterized once an initial point 
on the curve is chosen we conclude that S = AJF for some A > 0. This 
proves Kepler's first law. A comparison with the picture of the previous 
section shows that S = XT with A = —L/H. Indeed £ has foci and 
— AK/mL = K./mH = t and the long axis is equal to Xk/L = —k/H = 2a. 

It is not clear to us if Feynman was aware of his use of the existence 
and uniqueness theorem for differential equations. On page 164 of [S] he 
writes: "Therefore, the solution to the problem is an ellipse - or the other 
way around, really, is what I proved: that the ellipse is a possible solution to 
the problem. And it is this solution. So the orbits are ellipses." 

Apparently Feynman had trouble following Newton's proof of Kepler's 
first law. On page 111 of [3] the authors write: "In Feynman's lecture, this is 
the point at which he finds himself unable to follow Newton's line of argument 
any further, and so sets out to invent one of his own" . 

4 Newton's proof of Kepler's first law. 

In this section we discuss the original proof by Newton of Kepler's first law 
as given in [8|. The proof starts with a nice general result. 
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Theorem. Let S be a smooth closed curve bounding a convex region con- 
taining two points c and d. Let r{t) traverse the curve £ counterclockwise 
in time t, such that the areal speed with respect to the point c is constant. 
Likewise let r(s) traverse the curve £ counterclockwise in time s, such that 
the areal speed with respect to the point d is equal to the same constant. 

Let C be the tangent line to £ at the point r, and let e be the intersection 
point of the line M., which is parallel to C through the point c, and the line 
through the points r and d. Then the ratio of the two accelerations is given 
by 

d'^r d^Y 

Proof. Using the chain rule we get 

dr dv dt 

ds dt ds 

(Pv d^r ,dt.2 dr dH 
d^^d^'^d^' ^dt'd^' 
Because d'^r/dt^ is proportional to c — r and likewise d'^r/ds'^ is proportional 



|3 I II J|2 

e : r — c ■ r — d . 
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to d — r we see that 

(i^r iC^^r, fdt 2 ,d'^i^ dr (Pt dt 2< I'^^'^i /dt 2 

• l-^l = '^d^^di'd^'^d^^ ' ■ 'rft^l = (^) -Ir-el : |r-c| . 

Since the curve S is traversed with equal areal speed relative to the two 
points c and d we get 

,dr ,dr 
— • r — e = — ■ r — d 

dt ds 

and therefore also 

dt , , , 

— = r — e : r — d . 

ds 



In turn this implies that 

\2 I II II |3 I II j|2 

^ T — e : r — c = r — e : r — c • r — d 



which proves the theorem. □ 

We shall apply this theorem in case S is an ellipse with center c and focus 
d. Assume that r(t) traverses the ellipse S in harmonic motion, say 

d^r 
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Let b be the other focus of 8, and let f be the intersection point of the 
hne A/", passing through b and parallel to £, with the line through the points 
d and r. Then we find 

|d — e| = |e — f I , |f — r| = |b — r| 

which in turn implies that |e — r| is equal to the half long axis a of the ellipse 
8. We conclude from the formula in the above theorem that the motion in 
time s along an ellipse with constant areal speed with respect to a focus 
is only possible in an attractive inverse square force field. The converse 
statement that an inverse square force field (for negative energy H) indeed 
yields ellipses as orbits follows from existence and uniqueness theorems for 
solutions of Newton's equation F = ma and the above reasoning. This is 
Newton's line of argument for proving Kepler's first law. 

5 Conclusion. 

There exist other proofs of Kepler's law of ellipses from a higher view point. 
One such proof by Arnold uses complex analysis, and is somewhat reminis- 
cent to Newton's proof in Section 4 by comparing harmonic motion with 
motion under an force field [T]. Another proof by Moser is also very 
elegant, and uses the language of symplectic geometry and canonical trans- 
formations [6], [7]. However our goal here has been to give a proof which 
is as basic as possible, and at the same time is well motivated in terms of 
Euclidean geometry. 

It is hard to exaggerate the importance of the role of the Principia Math- 
ematica in the history of science. The year 1687 marks the birth of both 
modern mathematical analysis and modern theoretical physics. As such the 
derivation of the Kepler laws from Newton's law of motion and law of uni- 
versal gravitation is a rewarding subject to teach to freshmen students. In 
fact we were motivated for our work, because we plan to teach this mate- 
rial to high school students in their final grade. Of course, the high school 
students first need to get acquainted with the basics of vector geometry and 
vector calculus. But once this is achieved there is nothing in the way of 
understanding our proof of Kepler's law of ellipses. 

For physics or mathematics freshmen students in the university who are 
already familiar with vector calculus our proof in Section 2 is fairly short and 
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geometrically well motivated. In our opinion of all proofs it qualifies best to 
be discussed in a freshman course on mechanics. 
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